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ABSTRACT 

We present the results of high-resolution (1-0.4 A) optical spectroscopy of a sample 
of very low-mass stars. These data are used to examine the kinematics of the stars 
at the bottom of the hydrogen-burning main sequence. No evidence is found for a 
significant difference between the kinematics of the stars in our sample with I-K > 3.5 
(Msoi ft 12.8) and those of more massive M-dwarfs (Msoi ~ 7 — 10). A spectral atlas 
at high (0.4 A) resolution for M8-M9+ stars is provided, and the equivalent widths of 
Csl, Rbl and Ha lines present in our spectra are examined. We analyse our data to 
search for the presence of rapid rotation, and find that the brown dwarf LP 944-20 is 
a member of the class of "inactive, rapid rotators" . Such objects seem to be common 
at and below the hydrogen burning main sequence. It seems that in low-mass/low- 
temperature dwarf objects either the mechanism which heats the chromosphere, or 
the mechanism which generates magnetic fields, is greatly suppressed. 

Key words: stars: activity - stars: chromospheres - stars: kinematics - stars: low- 
mass, brown dwarfs - stars: rotation - techniques: radial velocities 



1 INTRODUCTION 

There have been suggestions in the recent literature of the 
existence of a young (Ss 10 8 yr) population of brown dwarfs 
in the solar neighbourhood masquerading as VLM stars. 
Hawkins & Bessell (1988) carried out a survey using UK 
Schmidt telescope plates, and claim that the reddest stars 
show a smaller (ie. younger) velocity dispersion than their 
mid- M- type stars. Kirkpatrick & McCarthy (1994) have 
shown that a linear extrapolation of their spectral-type ver- 
sus mass relation implies sub-stellar masses for spectral 
types later than M7. In such a situation, there would be a 
strong bias towards detecting only young, higher-luminosity 
(Msoi 13) VLM dwarfs. More compelling, however, are 
the results of the CCD survey of Kirkpatrick et al. (1994) 
which finds a significant excess of Mi > 12 (Mb i <; 12) 
dwarfs in the southern Galactic hemisphere: seven in 7.9 
square degrees versus none within the 19.4 area surveyed at 
northern latitudes. Given that the Sun lies ~ 30pc above 
the Galactic mid-plane, this would be consistent with the 
existence of a population of very low scale-height M-dwarfs. 
Kirkpatrick et al. suggest that this population is young ei- 
ther because such low-mass stars can only form at high 



Based on observations made at the European Southern Ob- 
servatory 3.6-m telescope, La Silla, Chile, and the 5-m Hale tele- 
scope, Palomar Observatory, USA. 



metallicities, or because the objects are substellar, and are 
detected only when young. These hypotheses can be tested 
by estimating a kinematic age for local VLM dwarfs, and our 
current analysis centres on applying that test. This work su- 
percedes the preliminary results presented by Reid, Tinney 
& Mould (1994 - hereafter RTM). 



2 OBSERVATIONS & ANALYSIS 

While proper motions are straight-forward (if tedious) to 
acquire for even very faint M-dwarfs, radial velocities are 
more difficult to measure. The low intrinsic luminosities of 
late-type dwarfs lead to only a relatively small number of 
those stars being accessible to accurate velocity measure- 
ment, particularly since those velocities should be obtained 
by cross-correlation on the stellar pseudo-continuum. The 
oft-used Ha emission line is produced via chromospheric 
activity and may not be a reliable velocity indicator. The 
present sample includes 37 stars or systems bright enough 
m the 7000-9000A region (1^17) to obtain radial velocities 
- 20 of these have M Bo i >12.8. 



2.1 ESO 3.6-m 

Spectra were obtained using the Cassegrain Echelle Spectro- 
graph (CASPEC) on the ESO 3.6m telescope on the nights 
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Table 1. Radial Velocity Measurements - Photometrically Selected Objects 
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NOTES: 

(a) Incorrectly referred to as BRI 0246-1703 by Tinney, Mould & Reid 1993. (b) These objects, though originally catalogued as 
proper-motion objects, were independently selected as late stars by an optical photometric survey (Kirkpatrick, Henry Sz Irwin 1997). 
They are therefore included in the "Photometric" sample, (c) A possible velocity variable, although the velocities are poor, 
(d) TVLM 513-42404A has a common proper motion companion with magnitude 1=18.56, I-K=4.35 (Tinney et al. 1995). The 
measurement of -9.5km/s was not used in calculating the mean, (e) Classified as a giant (M7-8III) by Kirkpatrick, Henry & Irwin 1997. 
Incorrectly referred to as BRI 2339-0447 by Tinney, Mould & Reid 1993. (f) I-K photometry is not available. A spectral type of M5.5 
has been assigned by Reid, Hawley & Gizis 1995. (g) G1473 is a known binary with separation 0.2-0.8" and P«16yr (Henry & 
McCarthy 1990). The components were not resolved by the spectrograph slit, (h) VB8 was used as a velocity reference at ESO, so ESO 
observations arc not included in the mean, (i) VB 10 was used as a velocity reference at Palomar, so Palomar observations are not 
included in the mean, (j) Optical astrometry indicates that G1831 may be an unresolved binary with a period of under two years. 
(Henry & McCarthy 1993 and references therein), (k) I— K photometry is not available. The spectral type is due to Kirkpatrick et al. 
1991. In view of the type mismatch between the object and template, this velocity should be treated with caution. (1) GL866AB is a 
known triple system, with two of the components having ~ 0.5" separation and P=2.20yr (Leinert et al. 1990). This probably explains 
the observed velocity variations. The components were not resolved by the spectrograph slit. The "mean" velocity shown should be 
interpreted with care in light of this. 
PHOTOMETRY REFERENCES: 

(1) Leggett 1992; (2) Tinney 1996; (3) Tinney, Mould & Reid 1993; (4) Kirkpatrick, Henry & Irwin 1997; (6) Hawkins & Bessell 1988; 
(7) Gullixson et al. 1995 and Harrington et al. 1993; (8) Weis 1988. 
POSITION REFERENCES: 

(1) Tinney 1993; (2) Tinney, Mould & Reid 1993; (3) Kirkpatrick, Henry & Irwin 1997; (4) Hawkins & Bessell 1988; (5) Gliese & 
Jahreiss 1991; (6) Luyten 1979; (7) Tinney 1996. 



of 1993 July 28-29 and 1994 February 8-9 (UT). The full 
width at half maximum (FWHM) resolution was 16km/s 
(eg. 0.43 A at 8000A), at a dispersion of 9 km/s per pixel. The 
wavelength range covered was 6400A to 9100A, though red- 
wards of 8045A the wavelength coverage is not complete and 
suffers inter-order gaps. For the faintest stars observed at 
ESO (1^16.3) exposures of up to 4 hours were used. Despite 



the fact that the signal-to-noise in an individual pixel at this 
resolution is low (S/Nw2.5 at 6500A, S/N«30 at 9200A), the 
long wavelength coverage of the instrument means that the 
final velocity precision is limited by the flexure, rather than 
S/N. Flexure is severe in CASPEC, so considerable care was 
taken in the observing and data reduction. Arc spectra were 
taken before and after each exposure, in order to provide 
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Table 2. Radial Velocity Measurements - Proper-motion 
Object V hel UT 
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NOTES, PHOTOMETRY and POSITION REFERENCES: See Table [[ 



a mean wavelength zero-point. The instrumental velocity 
system so constructed was checked by cross-correlating the 
wavelength calibrated night-sky spectra against each other - 
from which we derive a systematic uncertainty due to flexure 
of 2 km/s. 

The radial velocity system was calibrated initially us- 
ing observations of the preliminary IAU M-giant standard 
aCet {Vhel = -25.8 ±0.5 km/s; Trans.IAU,XXIB,275). By 
cross-correlation with a Cet we determined the heliocentric 
radial velocity of VB8 to be 14.5+1.0 km/s and of VB 10 
to be 34.4+1.3 km/s. VB8 was then adopted as the cross- 
correlation template for all of the CASPEC data, with an 
assumed Vhei = 14.5 km/s. Cross-correlations were carried 



out on an order-by-order basis, using regions of the spectrum 
free of atmospheric absorption and strong sky lines. The re- 
sults for each order were then combined as a weighted mean 
to produce a velocity shift for each exposure. Tables |l| and 
^| show measured radial velocities for each observation. The 
uncertainties quoted for individual observations are based on 
the uncertainty produced by the cross-correlation procedure 
of Tonry & Davis (1979)f|, combined with the systematic 

T RTM used an empirical technique to estimate cross-correlation 
velocity uncertainties, together with the expressions of Tonry & 
Davis (1979). Based on this it was stated that the Tonry & Davis 
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2km/s uncertainty due to spectrograph flexure - for most 
targets it is this systematic uncertainty which dominates. 



2.2 Palomar 5-m 

Radial velocities were obtained using the red arm of the Dou- 
ble Spectrograph on the Hale 5m telescope, on 1992 June 
19-22 (UT) and 1992 October 16-17 (UT). A 1200 line/mm 
grating was used in the following configurations: 76km/s 
(2.15A) resolution at a dispersion of 28km/s (0.80A) per 
pixel, over 8110-8735A in June 1992; and 74km/s (1.78A) 
resolution at a dispersion of 36km/s (0.8lA) per pixel over 
7906-7549 A. These are the same observations reported by 
RTM. However, in the course of checking our CASPEC ra- 
dial velocities against both the literature and the results of 
RTM, it became clear that there were problems with the 
RTM velocities, which at that time could not be checked 
against other observations owing to a lack of published ve- 
locities for late M-dwarfs. 

We have re-reduced the RTM data, taking the care to 
ensure a rigid velocity system with the CASPEC data. Cal- 
ibration arcs show flexure in the Double Spectrograph of up 
to 20km/s over the course of a night, so arcs were taken 
before and after each exposure of longer then 15 minutes. 
Cross correlations of arcs before and after long exposures 
reveal no flexure larger than 5km/s. We adopt a lower limit 
to the precision of our velocities of 3 km/s based on our abil- 
ity to calibrate out this flexure. 

VB 10 was adopted as the velocity template, with an 
assumed Vhel = 35 km/s (based on our CASPEC observa- 
tions). VB 10 was adopted in preference to VB8, as it was 
felt that at the lower spectral resolution of the Palomar data 
a good spectral-type match between template and object is 
more important. The results of this re-reduction are also 
shown in Tables |l] and |^. 



2.3 Combined Results 

The final velocities are weighted means of the individual 
observations. Amongst the current sample, only G1866AB 
shows compelling evidence for velocity variability. The ve- 
locity difference seen in exposures separated by only 2 days 
is 15.0±5.7 km/s. G1866AB has been resolved by infrared 
speckle techniques as a binary system with period of 2.2 
years and separation of between 0.2-0.8" by Leinert et al. 
(1990). However, Leinert et al. (1997) report a private com- 
munication from A.Duquennoy & D.Latham to the effect 
that a third component has been detected spectroscopically. 
Reid & Gizis (1997) also report the detection of an SB2 com- 
ponent, with a period inconsistent with that of the resolved 
components. Our current results add further support for the 
presence of a third component in this system. 

TVLM 832-10443 and TVLM 513-42404A both show 
possible evidence for velocity variability, although the data 
are of poor quality. In the case of TVLM 513-42404A where 



algorithm gave systematically low uncertainties. It is now known 
that this was due to an error in our implementation of their al- 
gorithm - the correct algorithm gives similar uncertainties to our 
empirical technique. 



Table 3. Radial velocities compared with other studies. 

Object Vhel Vhel AV hel 

Literature 



Vhel 
This work 



Marcy & Benitz (1989) 

G1699 -110.9±0.2 -112.3±5.1 1.4±S.l 

G1831 -57.1±1.0 -56.5±2.0 -0.6±2.2 

Delfosse, Forveille, Perrier & Mayor (1998) 

GJ1002 -41±1 -35.7±5.5 -5.3±5.1 

GJ1111 9±1 47.3±3.2 -38±3.3 

G1406 a 19.2±0.2 20.0±3.3 -0.8±3.3 

G1699 -lllil -112.3±5.1 -1.3±5.2 

Basri & Marcy (1995) 

GJ1111 11.4±3 b 47.3±3.2 -36±4 

BRI0021-0214 15.6±3 b 4.3±2.2 11.3±3.7 



NOTES: 

a - G1406 velocity is not that from Table 2 of Delfosse et al., 
but the more precise value quoted in Martin et al. 1997. b - 
Basri & Marcy used Gl 406 as a velocity standard with an 
assumed velocity of 16.5 km/s. The velocities they give have 
been corrected to using more precise velocity for Gl 406 given by 
Delfosse et al. 1998. 



three observations are available, we have rejected the dis- 
crepant value to estimate a mean velocity. For TVLM 832- 
10443 we use the mean of the two available velocities. 

In Table [] we compare the velocities derived in this 
study with those obtained by other studies - in particular 
those of Marcy & Benitz (1989), Delfosse et al. (1998) and 
Basri & Marcy (1995). An examination of this table shows 
that with only two exceptions our velocities are all consis- 
tent within 1-cr uncertainties with those obtained from the 
literature. 

Those exceptions are GJ 1111 and BRI 0021-0214. In 
the case of GJ 1111 we can find no plausible explanation for 
the velocity difference observed. The remote possibility that 
that GJ 1111 is a radial velocity variable is argued against by 
the consistent velocities for it determined at several epochs 
over the last 2 years at Observatoire de Haute Provence 
(X. Delfosse, priv.comm.). The target observed and the data 
reduction for this object have been throughly checked. We 
can find no fault to which this discrepancy (which corre- 
sponds to approximately one pixel on our detector) can be 
adduced. The observations of LHS2065, VB 8 and G1643 
made on the same night as GJ 1111, show that the velocity 
system for this night was consistent with that obtained on 
all other observing runs. 

Our velocities for BRI 0021-0214 show marginal evi- 
dence for variability. Moreover, they are also significantly 
different from that obtained by Basri & Marcy in November 
1993. Multiple observations on that run ruled out velocity 
variability on 90 minute timescales at the 3 km/s level (in- 
dicating the "system" is not a close, tidally-locked binary). 
The two combined sets of data therefore suggest ~ 2 — 12 
month velocity variability at the 10 km/s level. Further ob- 
servations are required to confirm or reject this hypothesis. 
It should also be noted that both we (§6) and Basri & Marcy 
find BRI 0021-0214 to be a rapid rotator (vsini « 40 km/s. 
This means its cross-correlation peak is significantly broad- 
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Table 4. Derived Velocities 



Object 


Vtan 






W a 




Sample c 


Astrometry 




km/s 


km/s 


km/s 


km/s 


km/s 




References 


BRI 0021-0214 


8.5± 0.5 


-8.8± 0.5 


-3.2± 1.5 


-2.6± 2.6 


9.7± 3.0 


A 


1 


TVLM 831-161058 


63. 6± 8.2 


-65.5± 6.0 


-43.2± 4.6 


10.6± 5.9 


79.1± 9.6 


A 


1 


TVLM 832-10443 


26. 3± 2.4 


9.5± 3.4 


-5.4± 2.8 


-23.9± 3.6 


26.3± 5.7 


A 


1 


LP 771-21/BR 0246-1703 


21. 8± 2.2 


16. 1± 3.0 


-24.0± 2.9 


12.3± 3.2 


31.4± 5.3 


A 


2 


LP 944-20/BRI 0337-3535 


10.3± 0.3 


-20.6± 0.8 


-15.9± 0.8 


-7.3± 0.7 


27.1± 1.3 


A 


2 


BRI 1222-1222 


26. 0± 1.9 


-20.8± 1.4 


-30.8± 1.5 


-18.2± 1.8 


41. 3± 2.8 


A 


2 


TVLM 513-46546 


4.0± 0.4 


-4.0± 1.3 


-12.5± 1.2 


1.5± 1.4 


13.2± 2.2 


A 


1 


TVLM513-42404A 


17.2± 5.1 


-30.6± 3.9 


-37.6± 3.4 


-34.2± 4.0 


59.3± 6.8 




1 


TVLM 513-8238 


24.2± 4.8 


-19.6± 3.4 


-36.2± 3.0 


-6.5± 3.5 


41. 7± 6.0 


A 


1 


TVLM 868-110639 


31. 4± 2.6 


-60.3± 3.7 


-28.7± 2.8 


-22.8± 3.6 


70.5± 5.9 


A 


1 


TVLM 890-60235 


17.7± 2.0 


7.9± 2.1 


-8.8± 2.1 


-0.5± 3.1 


11. 8± 4.3 


A 


1 


HB 2124-4228 


72.2±15.5 


-17.7± 8.1 


-82.8±10.0 


-6.5± 9.1 


84.9±15.8 


A 


1 


LHS 2/GJ1002 


45. 5± 0.8 


27.9± 0.7 


-49.6± 2.8 


16.1± 4.8 


59.2± 5.6 




2 


LHS112 


81. 5± 4.0 


-52.0± 2.7 


-50.1± 3.1 


-64.8± 3.5 


97.0± 5.7 




4 


ESO 207-61 


34.3± 2.9 


-61. 0± 2.3 


-118. 1± 2.0 


-29.4± 1.3 


136.1± 3.3 


B 


2 


LHS 234/G1 283B 


53. 0± 2.5 


45. 7± 1.8 


-10.6± 1.6 


21. 6± 1.9 


51. 7± 3.2 




4 


LHS 248/GJ1111 


22. 2± 0.4 


-57.0± 2.1 


-30.3± 1.7 


0.9± 1.7 


64.6± 3.2 




4 


LHS 2065 


22. 0± 0.3 


-23.1± 1.2 


-22.1± 0.9 


-21. 6± 1.1 


38.6± 1.8 


B 


3 


LHS 36/G1406 


53. 1± 0.3 


-36.9± 0.9 


-58.8± 1.6 


-18.9± 2.7 


71. 9± 3.3 




4 


LHS 333/G1473AB 


37.4± 0.7 


-43.0± 0.7 


-26.4± 0.9 


-6.2± 1.5 


50.9± 1.9 




4 


LHS 2351 


49. 0± 4.3 


-54.8± 2.9 


6.2± 2.9 


-2.0± 2.4 


55.2± 4.8 


B 


4 


LHS 2397a 


38. 3± 2.5 


-39.0± 2.1 


-51. 7± 1.3 


2.6± 1.4 


64.8± 2.9 


B 


3 


LHS 2875 


102.6±12.6 


-116. 5± 9.0 


-41. 9± 6.3 


2.3± 6.3 


123.9±12.7 




2 


LHS 2876 


102.6±12.6 


-118. 2± 9.0 


-41. 3± 6.3 


-0.2± 6.4 


125.2±12.8 


B 


2 


LHS 2924 


40.0± 0.6 


-1.8± 1.5 


-57.5± 1.7 


-28.4± 1.9 


64.2± 3.0 


B 


3 


LHS 2930 


37.1± 0.5 


-37.0± 1.6 


-37.5± 2.3 


-4.4± 1.8 


52.9± 3.3 


B 


3 


LHS 3003 


28. 8± 1.0 


-12.6± 1.3 


-37.1± 1.3 


-17.7± 1.3 


43.0± 2.2 


B 


2 


LHS 427/G1 643 


32. 7± 1.3 


11. 6± 1.3 


-39.8± 1.0 


4.8± 0.9 


41. 7± 1.9 




4 


LHS428/G1 644AB 


36. 2± 0.6 


13. 9± 1.4 


-42.9± 0.8 


6.0± 0.6 


45.5± 1.7 




4 


LHS429/G1 644C/VB8 


36. 0± 0.6 


9.4± 1.7 


-43.6± 1.0 


4.4± 0.7 


44.8± 2.1 




3 


LHS57/G1 699 


89. 4± 0.2 


-151. 3± 4.4 


-7.2± 2.3 


11. 9± 1.0 


152.0± 5.1 




4 


LHS 474/G1 752B/VB 10 


42. 4± 1.8 


44.9± 1.8 


-20.3± 1.5 


-11. 8± 1.1 


50.6± 2.7 


B 


1,2,3 


LHS511/G1 831 


45. 1± 1.8 


-73.3± 1.6 


-43.6± 1.3 


-0.7± 1.7 


85.3± 2.7 




4 


LHS 515 


267.9±55.0 


-107.2±37.2 


-126.1±29.8 


-234.9±26.7 


287.4±55.5 




4 


LHS 523 


55. 9± 3.4 


22.8± 2.2 


-57.5± 3.1 


-6.6± 3.5 


62.2± 5.2 




4 


LHS 68/G1 866AB 


52. 6± 0.6 


-73.9± 1.0 


-8.1± 1.5 


25.8± 2.2 


78.7± 2.9 




1 



NOTES: 

a - Note that the uncertainty ellipse associated with the measured space velocities does not have axes orthogonal to the (U,V,W) 
co-ordinate system. This means that the uncertainties in the U,V,W velocities are not independent, and so they do not add in 
quadrature to give the uncertainty in the total space velocity. 

b — Magnitude of the total space velocity, corrected for the basic solar motion as described in the text. 

c — Sample A is the set of photometrically selected objects with I-K>3.5 (ie Ms j>12.8). It also includes the binary companion 
TVLM 513-42404B, which is assumed to have the same radial velocity, distance and proper motion as TVLM 513-42404A. Sample B is 
the set of proper motion selected objects with I-K>3.5 (ie Mg ;>12.8). 
ASTROMETRY REFERENCES: 

(1) Tinney et al. 1995; (2) Tinney 1996; (3) Monet et al. 1992; (4) van Altena, Lee & Hoffleit 1995. 



ened, making it possible that the velocity variability ob- 
served is simply due to the uncertainty in measuring the 
centroid of a broadened peak. 

With the exception of these two objects, our velocities 
are consistent with external observations and we conclude 
that they are accurate to the stated precisions. In particu- 
lar, we are now confident that our re-reduction of the Palo- 
mar radial velocity data provides a robust velocity system, 
consistent with our ESO data and with independent obser- 
vations. The data presented in Tables [1] and therefore, 
supersede that presented in RTM, which should no longer be 
used. 



3 RADIAL VELOCITIES & SPACE 
VELOCITIES 

The current sample includes 37 stars; 20 fainter than Mb„i = 
12.8 (I-K > 3.5) and 17 brighter (the two stars with no I- 
K colour data are clearly more luminous than Mb i = 12.8, 
given their a spectral types of M5.5V and sdM2). A subset of 
these stars were identified originally through their substan- 
tial proper motions, and one might expect that subset to be 
biased toward higher velocities. To take this into account, 
we define two sub-samples: 

Sample A - is the set of photometrically selected ob- 
jects with I-K>3.5 (MboI > 12.8), containing 12 objects 
(Table [j]); 
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Table 5. Mean Velocities and Dispersions for VLM Stars. 



Weighted 

Sample A Sample B SampleA+B 
(km/s) (km/s) (km/s) 



Unweighted Early M dwarfs 

Sample A Sample B Sample A+B Hawley et al." 

(km/s) (km/s) (km/s) (km/s) 



<U> 
<V> 

<w> 



-13+7 
-20±6 
-8±4 



-34±14 
-42±11 
-13+4 



-17±8 
-30±6 
-10±3 



-18±7 
-27±6 
-8+4 



-34±14 
-42±11 
-13±4 



-25±7 
-33±6 
-10±3 



-18.6 
-32.1 
-13.6 



ajr 
a v 



18 
16 

23 



34 
30 
30 



26 
25 
25 



25 
22 
25 



45 
31 
34 



35 
28 
27 



38.2 
25.7 
20.5 



NOTES: 

a - The "Complete" sample results from Table 9 of HGR, but corrected for basic solar motion as for our data (§3). 



Sample B - is the set of proper motion selected objects 
with I-K>3.5 (Msoi > 12.8), containing 9 objects(Table ||). 



Table y lists the Galactic (U,V,W) space motions (U positive 
toward the Galactic Centre), where we correct for a Basic 
Solar Motion of (9,11,6) km/s (Mihalas & Binney 1981). In 
each case, the radial velocity and proper motion uncertain- 
ties were propagated into the (U,V,W) system. However, it 
should be noted that since the uncertainties in U,V and W 
are not independent, they do not add in quadrature to give 
the uncertainty in the total space velocity. 

None of the "confirmed" radial velocity variables in Ta- 
bles [l] and | (eg. G1866AB) are included in Samples A or B. 
The possible variables BRI 0021-0214 and TVLM 832-10443 
are included pending a final determination of their status 
- in any case their velocity variation is of the order of the 
quoted uncertainty in their velocities. The total number of 
objects in Samples A and B is only slightly larger than that 
used in the equivalent samples of RTM (ie 21 versus 18), 
but the precision is improved significantly and none of the 
current radial velocities are derived from Ha emission. Al- 
together, the conclusions we derive from these new data can 
be expected to be considerably more robust. 

Table |H] shows the means (<U>,<V>,<W>) and stan- 
dard deviations (au ,ov ,<jw) of the Galactic velocity distri- 
butions for Samples A, B and A+B. Both weighted and un- 
weighted estimates of these quantities are shown. The same 
quantities are shown for the 346 stars of the more luminous 
(My « 8 — 13, or Msoi £s 7.4 — 10.4) volume-complete sam- 
ple of Hawley, Gizis & Reid (1997 - hereafter HGR), with 
the same correction applied for basic solar motion. Figures 
[l]-^ show the Vtotai, U, V, and W velocity distributions for 
the same three samples, together with the corresponding 
distributions from HGR. 

Inspection of Table ^ shows that whether weighting is 
used or not makes no significant difference to the observed 
means - which are in all cases consistent with those observed 
for the more massive stars of HGR. The value of the velocity 
dispersion depends slightly on whether weighting is used, 
with the weighted values being somewhat smaller. Again, 
however, there is no indication in any case that the velocity 
dispersion for VLM stars is smaller than that seen for early 
M-dwarfs. 



Table 6. K-S Test Results - the probabilities, for each velocity 
component, that the hypotheses described in the text are true. 







Hypothesis 






(a) 


00 


(c) 


(d) 


Velocity 










U 


33% 


78% 


67% 


53% 


V 


33% 


89% 


77% 


32% 


w 


96% 


17% 


12% 


91% 



Because such an "eyeball" comparison can be mislead- 
ing, we have performed Kolmogorov-Smirnov tests (see eg. 
Press et al. 1986) to test the following hypotheses: 

(a) that Samples A and B are drawn from the same statis- 
tical population; 

(b) that Sample A+B is drawn from the same statistical 
population as the HGR sample; 

(c) that Sample A is drawn from the same statistical pop- 
ulation as the HGR sample; 

(d) that Sample B is drawn from the same statistical pop- 
ulation as the HGR sample; 

The results of these tests are shown in Table ^ and provide 
no evidence that samples A and B are drawn from different 
parent populations. Figs [l|j3] suggest possible kinematic dif- 
ferences between the VLM-dwarf and HGR samples. How- 
ever, the K-S tests show that all of the samples are indistin- 
guishable statistically. 



4 A SPECTRAL ATLAS AT 0.43A 
RESOLUTION 

High resolution echelle spectrographs efficient enough to ob- 
serve faint targets have been used now for several years 
to study the lowest mass stars (eg. Basri & Marcy 1995; 
Schweitzer et al. 1996). Despite this the available atlases 
for these objects are all based on data taken at low resolu- 
tion (eg. Kirkpatrick et al. 1991, 18A; Turnshek et al. 1985, 
8-12A). We therefore present in Figures ^J-^j an atlas show- 
ing CASPEC spectra of three objects; the M7 dwarf VB 8 
(Kirkpatrick et al. 1991); the M9 dwarf LHS2065 (Kirk- 
patrick et al. 1991) and the >M9V brown dwarf LP 944-20 
(Kirkpatrick et al. 1997; Tinney 1998). 
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Figure 1. The V( ota j, U, V and W velocity distributions for Sample A are shown as shaded. A solid Gaussian shows the weighted mean 
and dispersion. The dotted histogram is the "Complete" sample of HGR, normalised to the same total counts as Sample A. All velocities 
have been corrected for the basic solar motion as described in the text. 
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Figure 2. As for Fig. [j], but for Sample B. 



(km/s) 



These spectra have been flux calibrated using obser- 
vations of the standard star fj, Columbae (Turnshek et al. 
1990). Because our observations were not made with a spec- 
trophotometric slit, we only obtain a relative inter-order 
calibration. An absolute calibration has been obtained by 
scaling to the published I-band photometry of these objects 
(Leggett 1992) using a zero-magnitude flux of 2550 Jy (Reid 
& Gilmore 1984). The overall flux calibration is accurate 
to ±10% - however the inter-order calibration is much bet- 
ter and good to a few percent. Photon counting errors were 
carried throughout the reduction. The signal-to-noise ratio 
(SNR) as a function of wavelength for our spectra is sum- 
marised in Fig [lo| Digital copies of the atlas spectra (includ- 
ing the SNR spectra) have been deposited with the NASA 
Astrophysics Data Centre. 



4.2 Molecular Features 

The dominant molecular features in these spectra are due 
to TiO and VO. The figure shows the locations of identified 
bandheads due to these two species. Also shown are the spec- 
troscopic "systems" to which these bandheads have been as- 
signed. Information on the actual state transitions to which 
these systems correspond can be found in the references pro- 
vided below - in particular see Pearse & Gaydon (1976). It 
is important to emphasise that these identifications repre- 
sent only the bandheads - more widely spread transitions 
cover almost all of the optical and infra-red spectra of these 
stars. Between 8700A and 8850A for example there are no 
major bandheads, but the spectra clearly contain numerous 
molecular features, as can be seen in the laboratory spectra 
of Gatterer et al. (1957). 



4.1 Terrestrial Features 

The spectra shown have not been corrected for terrestrial 
absorption due to O3 and H2O (cf. Allen 1976, §58 and 
references therein). The locations of terrestrial absorption 
in these spectra can be seen in the panel showing the nor- 
malised flux from the standard star /i Col. The effects of O3 
are quite regular and straightforward to disentangle from 
features in the stars themselves. This is not true for H2O ab- 
sorption which is complex and irregular. We therefore show 
in the stellar panels the locations of the strongest lines of 
H2O absorption. Also shown are the locations of the stellar 
absorptions present in fi Col itself. 



4.2.1 TiO 

All the TiO bandheads are observed to be degraded to the 
red. Where no reference is given consult Pearse & Gay- 
don (1976), Solf (1978), Gatterer et al. (1957) and Turn- 
shek et al. (1985). Indicated are: the 7 or "Red" system 
bandheads at 6651.5, 6681.1, 6714.4, 6719.3, 6746.7, 6781.3, 
6814.7, 6849.9, 6852.3, 6883.6, 6918.9, 7054.5, 7059.2, 
7087.9, 7093.1, 7124.9, 7125.6, 7130.4, 7159.0, 7197.7, 

7219.4, 7230.8, 7269.0, 7589.6, 7628.1, 7666.4, 7672.1, 
7705.2, 7742.8, 7751.6, 7782.7, 7820.1, 7828.0, 7861.0, 
7881.2, 7907.3, 7948.6, 8205.8, 8250.6, 8289.0, 8302.9, 

8334.5, 8375.5, 8386.5, 8419.5, 8442.3, 8451.7, 8457.1, 

8471.6, 8505.5, 8513.1, 8558.4, 8569.4, and 8645.5A; the 7/ 
or "Orange-Red" system bandheads at 6569.3, 6596.2 and 
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Figure 3. As for Fig. [l], but for Sample A+B. 



6629. OA; an unnamed system with bandheads at 8432, 8442 
and 8452A (Solf 1978); and an "Infrared" system with band- 
heads at 8859.6, 8868.5, 8937.4, 8949.8, 9014.6 and 9094.5A 
(Philips 1950). Also shown is a bandhead which fits into no 
system at 8198. 5A (Solf 1978) 



4.2.2 VO 

Unlike TiO which produces distinctive bandheads degraded 
to the red, VO in M-dwarfs produces more diffuse absorp- 
tion. Where no reference is given consult Pearse & Gaydon 
(1976), Solf (1978), Gatterer et al. (1957) and Turnshek et 
al. (1985). Indicated are: the B or "Red" system bandheads 
at 6876.2, 6894.0, 6919.0, 6951.6, 7011.0,7070.2, 7131.7, 

7418.2, and 8642. lA; the C or "Infrared" system bandheads 
at 7333.8, 7345.3, 7372.4, 7374.8, 7393.2, 7405.2, 7433.6, 

7435.3, 7444.0, 7454.0, 7466.0, 7472.1, 7492.0, 7534.1, 
7850.9, 7865.0, 7896.0, 7899.6, 7918.4, 7928.5, 7938.9, 

7947.7, 7960.1, 7967.2, 7973.1, 7982.1, 8520.9, 8537.7, 

8572.8, 8575.3, 8590.7, 8597.2, 8604.0, 8624.0, 8648.6, 8657.9 
and 8666.6A. 
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Figure 10. Signal-to-noise ratios as a function of wavelength for 
the spectra shown in Figs WS. 



KI weakens slightly with decreasing effective temperature, 
or that VB 10 has an unusually high metallicity or gravity. 



4.2.3 FeH 

FeH bandheads are seen degraded to the red. Indicated 
are the 4 A- 4 A system bandheads at 7786, 8692 and 9020A 
(Philips et al. 1987). 

4.3 Atomic Features 

Atomic features shown are: Ha 6562. 8A; Li I 6708.0A; KI 
7664.9A and 7699.0A; Rbl 7800. 2A and 7947.6A; Nal 
8183.3A and 8194.8A; Csl 8521.4A and 8943. 6A; Bal 
7911. 3A. With the exception of Ha all these lines are res- 
onance transitions. The weak lines of Call, Mgl and Til 
seen in earlier M-dwarfs (Kirkpatrick et al. 1991) are un- 
detectable amongst the morass of molecular absorptions in 
these M7-M9+ dwarfs. The Li I detection seen in LP 944-20 
in Fig. ^ is discussed in Tinney (1998). 

It is particularly interesting to note the pronounced 
width of the K I lines in VB 10 compared with the two cooler 
objects. Schweitzer et al. (1996) have modelled the atomic 
lines of VB10, and find that this extreme strength of KI 
could be due to either high metallicity ([M/H]=+0.5) or 
high gravity (logg=5.5). Schweitzer et al. did not explore 
the temperature dependence of KI line strength, however 
the spectra shown in Fig. would seem to imply that either 



5 SPECTRAL FEATURES 
5.1 Atomic Lines 

In this section we deal only with the high resolution obser- 
vations made at La Silla. Na I and K I lines are prominent in 
almost all of these spectra, however measurement of them 
is made difficult by the presence of atmospheric absorption, 
which we did correct for this study. 

Csl and Rbl lines, however, lie in regions of the spec- 
trum clear of atmospheric absorption, and become promi- 
nent for late-type dwarfs (Tinney 1997). The measured 
equivalent widths of these lines are shown in Table [?]. Be- 
cause no part of the spectrum is free of stellar molecular 
absorption, the assignment of a true continuum is impossi- 
ble. As a result whenever we refer to equivalent widths, they 
are always defined relative to the apparent continua avail- 
able near the lines. The uncertainties presented in the table 
are based on photon-counting errors and do not include sys- 
tematics due to pseudo-continuum placement. However, our 
measurements should be self-consistent within this set of ob- 
servations and with other observations obtained at similar 
resolution. We have also included results from Martin et al. 
1997 for two cool objects discovered in the DENIS survey 
(Delfosse et al. 1997; Tinney, Delfosse & Forveille 1997). 
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Figure 4. A spectral atlas for the objects VB8, LHS2065, and LP 944-20 at 0.43A resolution. The normalised spectrum of the standard 
star /iCol shows the location of terrestrial atmospheric features. The bottom of the scale in each panel is 0.0. 
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Figure 5. A spectral atlas for the objects VB8, LHS 2065, and LP 944-20 at 0.43A resolution (continued). 
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Figure 6. A spectral atlas for the objects VB8, LHS 2065, and LP 944-20 at 0.43A resolution (continued). 
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Figure 7. A spectral atlas for the objects VB8, LHS 2065, and LP 944-20 at 0.43A resolution (continued). 
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Figure 8. A spectral atlas for the objects VB8, LHS 2065, and LP 944-20 at 0.43A resolution (continued). 



© 1998 RAS, MNRAS 000, 



C.G. Tinney & I.N. Reid 




I S'O 002 00 I 0£ 02 1 0> Of 02 OL 



Figure 9. A spectral atlas for the objects VB8, LHS 2065, and LP 944-20 at 0.43A resolution (continued). 
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Table 7. Cs I and Rb I Observations 



Table 8. CASPEC Ha Emission Measurements 



Object 

LP 944-20 
BPJ 1222-1222 
ESO 207-61 
LHS 2065 
LHS 2351 
LHS 2397a 
LHS 2875 
LHS 2876 
LHS 3003 
LHS 427 
LHS 428 
LHS 429 
LHS 474 

DENIS-P J1058 6 
DENIS-P J 1228 6 



I-K 

4.58 
4.31 
4.18 
4.46 
3.65 
4.26 
2.05 
3.59 
3.66 
2.30 
2.16 
3.42 
3.97 

5.1 
5.5 



Pseudo-Equivalent Widths (A) 
RbI 7800A RbI 7948A Csl 852lA 



0.87±0.05 
0.88±0.03 
0.94±0.06 
0.81±0.03 
0.76±0.05 
0.87±0.02 
0.16±0.02 
1.20±0.20 
0.84±0.02 
0.23±0.01 
0.19±0.01 
1.02±0.02 
0.93±0.05 



0.73±0.05 
0.83±0.05 
1.03±0.14 
0.62±0.05 
0.63±0.05 
0.96±0.05 

1.20±0.10 
0.86±0.03 
0.17±0.01 

0.95±0.03 
0.93±0.05 



0.48±0.02 
0.39±0.03 

_a 

0.33±0.05 
0.30±0.03 
0.50±0.10 



0.34±0.02 

<0.03 

<0.03 

0.36±0.03 

0.37±0.07 

2.1 

4.3 



NOTES: 

a - Csl 852lA was detected in ESO 207-61, but could not be 
measured because its apparent radial velocity moved the 
CsI852lA line onto the end of an echelle order. 
b - I— K from Delfosse et al. 1997. Equivalent widths for Cs I from 
Martin et al. 1997. Mj, D ; estimates from Delfosse et al. 1997. 




Figure 11. Strength of Csl lines as a function of I— K colour for 
the objects in Table ^. 



Cs I line strengths are plotted in Figure [H| The ob- 
servations describe a plateau at EW=0.2-0.3A throughout 
the late M-dwarf range, with considerably weaker absorption 
for the early M-dwarfs (as represented by the two upper lim- 
its), and considerably stronger absorption for objects below 
~2000K. This diagram highlights the possible utility of the 
Csl 8521 A line as a low-resolution spectroscopic indicator 
of effective temperature - at least for the purpose of dis- 
criminating extremely cool objects like DENIS-P J1228 and 
DENIS-P J1058 (for which it has been proposed that an en- 
tirely new spectral type ("L") is required; Kirkpatrick 1998) 
from late-type M dwarfs. 



Object 

BRI 0021-0214 
LP 944-20 

BRI 1222-1222 
ESO 207-61 

LHS 2065 
LHS 2351 
LHS 2397a 

LHS 2875 
LHS 2876 
LHS 3003 
LHS 427 
LHS 428 

LHS 429 



LHS 474 



UT 



1993/6/28 
1993/6/28 
1994/2/09 
1994/2/08 
1994/2/08 
1994/2/09 
1994/2/08 
1994/2/09 
1994/2/08 
1994/2/09 
1994/2/09 
1994/2/09 
1994/2/09 
1994/2/08 
1993/6/28 
1994/2/08 
1993/6/28 
1994/2/08 
1994/2/09 
1993/6/28 



07:55:22 
09:53:58 
01:15:52 
09:00:17 
02:17:06 
03:04:46 
04:22:29 
05:44:56 
06:13:18 
06:59:29 
07:59:21 
08:40:30 
07:35:41 
09:37:38 
03:04:22 
09:42:51 
02:17:11 
09:25:37 
09:31:43 
03:32:50 



Ha EW 
(A) 

<2 
<2 
0.6±0.2: 
4.7±0.5 
2.8±0.5 
2.6±0.8: 

24.9±0.3 
6.1±0.1 

47.3±0.8 

34.6±1.0 
0.20±0.02a 
1.6±0.5 
4.1±0.3 
0.24±0.01a 
1.8±0.4 
1.6±0.1 

20.4±1.3 
3.7±0.2 
4.1±0.1 
4.4±0.7 



[ L Hn ]b 

L Lbol J 



-5.1 
-5.1 
-5.6 
-4.7 
-4.9 
-5.0 
-4.0 
-4.3 
-3.7 
-3.9 

-4.9 
-4.6 

-3.5 
-3.6 
-3.7 
-4.5 
-4.4 
-4.7 



NOTES: 

a - "<" indicates a 3-cr upper limit. ":" indicates a marginal 
detection, "a" indicates the line was detected in absorption. 
Other results are actual measurements. 

b - Upper limits are indicated by <. Measurements have an 
uncertainty of ±0.2 largely due to the calibration from 
equivalent width to [Ljici/L;, ;]. 



m o 
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log 10 (S c /F bol ) = 0.076 - 2.35(I-K) + 0.25(I-K) 2 



■* . -P. . 



3 
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Figure 12. Ha "bolometric correction" as a function of I— K 
colour for VLM stars from Tinney, Mould & Reid 1993. Also 
shown is the calibration relation adopted to derive the [Lu a /Li, ;] 
values derived in Table ^ 



5.2 Ha Emission 

Ha emission is seen in low mass stars when magnetic fields 
heating the outer atmosphere producing temperatures of up 
to ~10 6 K in the outer corona. Closer to the photosphere the 
chromosphere shows lower temperatures of 10 4 K, where 
resonance lines of Call and Mgll, and the hydrogen Balmer 
lines are produced (HGR). Ha emission equivalent widths 
for our Palomar observations were reported in RTM. Table 
^| lists similar measurements for our ESO program objects, 
including only those observations yielding either detections 
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or significant upper limits. Because of poor seeing during 
our 1993 July observations, these come mostly from 1994 
February. 

Equivalent width (EW) is a measure of line flux (Fhc), 
measured relative to the local continuum flux density (S c ). 

EW« 

Sc 

A large equivalent width does not necessarily imply strong 
chromospheric activity. A better indicator is the ratio of the 
luminosity in the Ha line to the total bolometric luminosity 
of the star, Lua/i^boi, or its common logarithm([LHa/Li, oi ]). 
Since our observations were not made with a spectrophoto- 
metry wide slit, line fluxes could not be measured directly. 
However, the EW can be converted to L/Ha/Lboi if the ratio 
of the local continuum flux density S c to the bolometric flux 
(Fboi) is known. Since, 

L Ha = Fh« = EW.S C 
Lboi Fboi Fboi 

The ratio S c /Fboi is just a bolometric correction. We 
can therefore derive a functional form for it, dependent on 
effective temperature or bolometric luminosity. We have cal- 
culated this quantity for the sample of objects with cali- 
brated optical spectra and bolometric fluxes given in Tinney, 
Mould & Reid (1993), and the results are shown in Fig Pl2l 
The mean calibration was then used to derive the [LHa/L^j] 
estimates shown in Table || with an estimated uncertainty 
of ~ ±0.2. Full discussion of these results is deferred to §7. 



6 ROTATION 

Estimates of the rotational broadening of our CASPEC tar- 
get stars have been obtained by cross-correlating the tar- 
get spectra against a template (cf. Tonry & Davis, 1979). 
The dispersion^] of the cross correlation peak (/i) is then the 
quadratic sum of the dispersion of the instrumental broad- 
ening in each spectrum (r - assumed the same in both), any 
intrinsic broadening in the template spectrum (ero) and the 
rotational broadening of the target spectrum (a). Uncertain- 
ties in the measure ment of [i render it impossible to mea- 
sure values of a < */ 2r 2 + <r§. In the case of our CASPEC 
observations, r = 6.8km/s, setting a detection limit for ro- 
tational broadening (usini) of 9.6 km/s. The Palomar data 
have a FWHM resolution of 75 km/s and set no useful limits 
on stellar rotation. 

Following Tonry & Davis we estimate the width w of 
the cross-correlation peak n by fitting a parabola. However, 
unlike Tonry & Davis we have not sought to calibrate the re- 
lationship between this parabola width w and the Gaussian 
dispersion fj,. Rather we derive an empirical calibration be- 
tween w and usini by cross-correlating "spun-up" spectra of 
VB 8 against an independent observation of the same star. 
The rotationally-broadened spectra are constructed using 
the profiles given by Gray (1992, equation 17.12, e = 0.5), 

t In this section we use the term "dispersion" to describe the 
quantity a in the usual Gaussian parameterisation f(x) oc 
e — {."/") ^ w hich is related to the full width at half maximum 
(FWHM) by <t FWHM/2.354. 



Table 9. CASPEC Rotational Measurements 



Object 






I IT 




vs\Tii a 












(km/s) 


BRI 0021-0214 


28 


Jul 


1993 


07:55:22 


42.5±8.4 




29 


Jul 


1993 


08:36:41 


40.1±8.6 


LP 944-20 /BRI 0337-3535 6 


28 


Jul 


1993 


09:53:58 


28.3±1.9 




29 


Jul 


1993 


10:17:28 


28.4±3.9 




09 


Feb 


1994 


01:15:52 


28.1±1.8 


BRI 1222-1222 


08 


Feb 


1994 


09:00:17 


<17 


BR 2339-0447 d 


28 


Jul 


1993 


06:22:05 


<22 


ESO 207-61 


08 


Feb 


1994 


02:17:06 


<18 




09 


Feb 


1994 


03:04:46 


<19 


LHS 234/G1 283B 


09 


Feb 


1994 


04:25:04 


<14 


LHS 2065 


08 


Feb 


1994 


04:22:29 


<18 


LHS 2351 


09 


Feb 


1994 


05:44:56 


<16 


LHS 2397a 


08 


Feb 


1994 


06:13:18 


23.1±3.2 il 




09 


Feb 


1994 


06:59:29 


23.6±1.8 6 


LHS 2875 


29 


Jul 


1993 


23:42:57 


<22 




09 


Feb 


1994 


07:59:21 


<17 


LHS 2876 


09 


Feb 


1994 


08:40:30 


<18 


LHS 3003 


09 


Feb 


1994 


07:35:41 


<13 


LHS 427/ Gl 643 


28 


Jul 


1993 


02:50:07 


<19 




08 


Feb 


1994 


09:37:38 


<18 


LHS428/G1 644AB 


28 


Jul 


1993 


03:04:22 


30.5±3.6 C 




08 


Feb 


1994 


09:42:51 


<19 


LHS 474/G1 752B /VB 10 


28 


Jul 


1993 


03:32:50 


<19 




29 


Jul 


1993 


04:25:38 


<19 


LHS511/G1 831 


28 


Jul 


1993 


05:47:10 


<19 



NOTES: 

a - < indicates a 3-ct upper limit calculated as described in the 

text. Other results are actual measurements. 

b - The velocities shown are only marginally above the 3-cr 

detection threshold for these observations. 

c — Evidence was seen in the cross-correlation function for a 

double peak. This system is multiple so the width measured is 

probably due to orbital motion, not rotation. 

d - An M7-8III giant (Kirkpatrick, Henry & Irwin 1997). 

computed for a series of velocities between and 50 km/s. 
The resulting w-to-vsini calibration was calculated indepen- 
dently for each order. 

Our calibration assumes that VB 8, the template for 
our VLM sample, has negligible rotation and there are no 
previous observational constraints of this parameter. As a 
test, we have repeated our w-usim calibration procedure 
using G1643, which has a measured rotational upper limit 
of usini < 2.7km/s (Delfosse et al, 1998), cross-correlating 
against the VB8 template. The resulting relation is essen- 
tially identical to that derived using VB 8 alone, implying 
that the latter star has a rotation substantially lower than 
our instrumental resolution of r = 6.8 km/s. 

The rotational velocity for each target star was then es- 
timated as follows: each order was cross-correlated with the 
appropriate template star order; the width w of the highest 
peak in the cross-correlation was measured, if its confidence 
level exceeded 90% (Tonry & Davis, equation 33); the uncer- 
tainty in the measurement of this width Aw was also esti- 
mated (Tonry & Davis, equation 25) ; if the measured width 
was less than ii)o + 3Ara, then this measurement is treated as 
an upper limit at the usim corresponding to w = wo + 3Aw, 
otherwise the measured width was treated as a velocity mea- 
surement at the appropriate usim; and lastly, the resulting 
rotation measures were averaged over the available orders. 
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The results of this process are shown in Table Typical 
uncertainties in Aw range from 0.2km/s for bright stars to 
2km/s for faint ones. This results in 3-cr rotation detection 
thresholds ranging from 15 to 22 km/s. In other words, given 
the resolution of our observations, we are only sensitive to 
very fast rotators. Significant rotation was detected in only 
two objects: BRI 0021-0214 and LP 944-20. LHS 2397a seems 
to show possible evidence for rotation right at the detection 
limit for this star of ~ 23 km/s. G1644AB showed some evi- 
dence for broadening in 1993 July, but none in 1994 Febru- 
ary. Given the multiple nature of this system, and evidence 
seen for a double peak in the cross correlation function, it 
seems reasonable to assume that the broadening observed 
in 1993 is more likely due to orbital motions than stellar ro- 
tation. The effect of rapid rotation can be seen in the spec- 
tra for LP 944-20 plotted in Figures ^||], where both atomic 
and molecular features are noticeably broader than those 
observed in LHS 2065 or VB 10. Of the two stars in which 
we detect rotation, only BRI 0021-0214 has been previously 
measured. The value of 41±8 km/s we measure for BRI 0021- 
0214 compares favourably with the 40±7km/s measured by 
Basri & Marcy (1995). 



7 ROTATION & ACTIVITY 

The standard picture of the chromospheric activity (as re- 
vealed by Ha emission) is that it is a product of stellar mag- 
netic fields, via the effects of stellar rotation and convection 
combining to form a dynamo. In this picture, the activity is 
related directly to stellar rotation. Such a model is now well 
established for F,G and K-type stars (Hartmann & Noyes 
1987). Stellar rotation rate is also a function of age - the 
rotation of a star slows as it loses angular momentum and 
"spins down". The spin down rate is generally held to be a 
function of mass, with lower mass stars spinning down more 
slowly (eg. Stauffer et al. 1994). 

As a result a general relationship between chromo- 
spheric activity and age can be deduced, at least for F, G 
and K dwarfs, though this relationship is complicated by 
the fact that stars are born with an initial spread in an- 
gular momentum, producing considerable scatter about the 
overall relationship between chromospheric activity and age. 
This is evidenced by the large spread seen in the rotational 
velocities of stars in young and intermediate-age clusters: for 
example, G dwarfs in the Pleiades (~100 Myr) have vsini <6 
- 50 km/s, while M dwarfs in the Hyades («600Myr) show 
vsini <10 - 20 km/s (Stauffer et al. 1994). 

The situation for low mass M dwarfs, however, is less 
clear. The appearance of Ha emission increases with spectral 
type from M0 (« 10%) to M6 (« 50%). It has been suggested 
that this is because Ha becomes easier to detect as the con- 
tinuum flux at 6562A decreases. However, the ratio of the Ha 
luminosity to the bolometric luminosity (LHa/iboi) spans 
an approximately constant range of activity at [LHa/Lboi] 
= —3.5 to —4.2 (HGR), rather than declining with decreas- 
ing luminosity, as one would expect if a selection effect were 
responsible for the higher dMe fraction. 

The straightforward interpretation of this result is that 
the proportion of rapidly rotating M-dwarfs increases with 
spectral type. Delfosse et al. (1998) have studied the rota- 
tion of a sample of M-dwarfs, and find the expected correla- 



tion with rotation - all stars with vsini > 15 km/s show 
chromospheric activity at the [Lh q / Lboi] = —3.5 to —4 
level, while stars with no significant rotation detected show 
logio [LHa/Lboi] = —4 to —4.5. Thus, this result suggests 
that the higher fraction of late-type dMe dwarfs may reflect 
longer spin-down times at very low masses. 

However, the fast rotators detected in our CASPEC 
study throw a spanner into this neat picture. Basri & Marcy 
(1995) first pointed out that BRI 0021-0214 is extremely late 
in spectral type (>M9.5V, Kirkpatrick et al. 1997), and a 
fast rotator, but shows no Ha in emission. Since then Tin- 
ney et al. (1997) have found BRI 0021-0214 to emit Ha at 
very weak levels (1.30±0.05A EW). This corresponds to only 
[LHa/Lboi] = —5.3 (cf. §5.2) making this to all intents and 
purposes a chromospherically inactive star, and directly con- 
tradicting the generally held activity-rotation connection. 

LP 944-20, with [L HQ /L b oi] ~ -5.6 and vsini = 
28 km/s, represents the detection of a second member of 
this class of "inactive, rapidly rotating" objects. Moreover, 
in this case we can add the further physical parameters that 
the inactive, rapid rotator is a brown dwarf of mass 0.06Mq , 
and age 475-650 Myr (Tinney 1998). The recently discov- 
ered brown dwarf DENIS-P J1228. 2-1547 (wsim=20 km/s, 
[LHa/Lboi]~-5.3), and the brown dwarf candidate DENIS- 
P 1058.7-1548 (23 km/s, «-5.4) also fall into this group 
(Martin et al. 1997; Tinney et al. 1997). These objects may 
represent an entire class of objects at or below the hydrogen 
burning main sequence which show substantially reduced 
levels of chromospheric activity, despite the presence of lev- 
els of stellar rotations which would produce significant ac- 
tivity in an MO-6 dwarf. 

Is this violation of the rotation-activity connection 
purely a function of mass? Observations of nearby clusters 
suggest that this may be the case. Activity (as measured by 
[LHa/Lboi]) amongst Pleiades M-dwarfs rises monotonically 
with decreasing luminosity (and mass) until Mtoi ~ 10, at 
which point the relative flux in Ha drops sharply. Figure 
6(b), for example, of Reid, Hawley & Mateo (1995) shows 
that for luminosities below Mtoi ~ 10, [LHa/Lboi] falls from 
~ —3.5 to ~ —4.5. Similar results have been derived by 
Stauffer et al. (1994), Hodgkin, Jameson & Steele (1995) 
and Oppenheimer et al. (1997), while Reid & Hawley (1998) 
have detected comparable behaviour amongst Hyades stars. 

Spectroscopy of Pleiades brown dwarfs shows that sim- 
ilar mass objects to LP944-20 and DENIS-P J 1228 (Rebolo, 
Zapatero Osorio, Martin 1995; Basri, Marcy & Graham 
1996; Zapatero Osorio et al. 1997) have Ha lines with equiv- 
alent widths of ~ 5 — 10A, which corresponds to [LHa/Lboi] 
» -4.5 to -4.8 (§5.2) at age ~ 70 - 120 Myr. 

We therefore conclude that (1) we see evidence in the 
Pleiades for a trend towards lower chromospheric activity 
at masses below « 0.1 Mq, and (2) that we see no evidence 
for a significant change in this trend as the brown dwarf 
limit is crossed. Moreover, the Pleiades VLM stars studied 
by Oppenheimer et al. (1997) are all rapid rotators, with 
vsini in the range 37-65 km/s. This implies that Pleiades age 
objects just above the H-burning limit already violate the 
rot at ion- activity connection. The "inactive rapid-rotators" 
detected in the field continue this violation below the H- 
burning limit. Age does not seem to play a determining role, 
since both ~ 70 - 100 Myr VLM stars and ~ 500 - 1000 Myr 
brown dwarfs are found to rotate rapidly, and yet to be 



© 1998 RAS, MNRAS 000, 0-|l| 



18 (7. 67. Tinney & I.N. Reid 



substantially inactive. However, there is evidence for a de- 
crease in activity as VLM stars and brown dwarfs age, since 
[L/Ha/Lboi] falls from ~ —4.5-— 4.8 to ~ —5.1 — 5.3 as we 
move from Pleiades age brown dwarfs to ~ 500 — 1000 Myr 
brown dwarfs like LP 944-20 and DENIS-P J1228. 2-1547. 

It should be noted however, that the possible identi- 
fication of decreasing mass as a parameter well correlated 
with the violation of the rotation-activity connection, does 
not necessarily imply that mass is the controlling physical 
parameter. The absence of chromospheric activity implies 
either the absence of magnetic fields, or the suppression of 
the chromospheric heating mechanism. In either case this 
could be physically due to luminosity or effective temper- 
ature dependent effects (eg. the formation of dust which 
begins at VLM effective temperatures), as well as to the 
internal structure imposed on these objects by their mass. 

Confirmation of the rapid rotation of these objects via 
the search for periodic modulation needs to be made a high 
priority. Many very low-mass stars show periodic modula- 
tions due to star spots. Even if the inactive, rapid rotators 
have no spots, they may exhibit periodic variation due to 
the transit across their disk of dust clouds. In either case 
(spots or clouds) the result will be the passage of regions of 
apparently low effective temperature across their disk. Such 
regions should show a clear signature in the regions of the 
spectrum where TiO absorption is strong in late M-dwarfs, 
but absent in cooler objects like DENIS-P J1228 (Tinney et 
al. 1997). 

The inactive, rapid rotators represent a direct counter- 
example to the generally held picture of the activity-rotation 
connection. Unravelling of this mystery will provide valuable 
insights into both the interiors of brown dwarfs, and the 
mechanics of magnetic field and chromosphere generation. 



8 CONCLUSION 

We have obtained high resolution optical spectroscopy of a 
sample of VLM stars and one brown dwarf. Analysis of the 
kinematics of these stars show no evidence for a significant 
kinematic difference between either our photometrically- 
selected or proper-motion-selected sample, and a sample of 
more massive nearby M-dwarfs. 

We present 0.43A resolution observations of two VLM 
stars and one intermediate age brown dwarf, for use 
as a spectroscopic atlas, together with identifications for 
their major stellar atomic and molecular features. Pseudo- 
equivalent widths for the Csl and Rbl features present 
in these are obtained. The Csl 8521 A feature in particular 
shows promise for discriminating L-type from M-type dwarfs 
in low resolution spectra. 

Measurements of the Ha emission equivalent widths 
and stellar rotations are presented. These reveal that the 
brown dwarf LP 944-20 is an "inactive, rapid rotator" , like 
the M-dwarf BRI 0021-0214 and the VLM Pleiades objects 
of Oppenheimer et al. (1997). This discovery shows that vi- 
olation of the rot at ion- activity connection continues below 
the bottom of the H-burning main sequence. The fact that 
such violation occurs for both Pleiades-age VLMs and for 
a ~500Myr old brown dwarf, indicates that age is not the 
determining factor in such a violation. 
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